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Nonexistence of simple hyperbolic blow-up 
for the quasi-geostrophic equation 

By Diego Cordoba 
In memory of Maria del Carmen Gazolaz 
1. Introduction 

The problem we are concerned with is whether singularities form in finite 
time in incompressible fluid flows. It is well known that the answer is "no" in 
the case of Euler and Navier-Stokes equations in dimension two. In dimension 
three it is still an open problem for these equations. 

In this paper we focus on a two-dimensional active scalar model for the 
3D Euler vorticity equation. Constantin, Majda and Tabak [7] suggested, 
by studying rigorous theorems and detailed numerical experiments, a gen- 
eral principle: "If the level set topology in the temperature field for the 2D 
quasi-geostrophic active scalar in the region of strong scalar gradients does not 
contain a hyperbolic saddle, then no finite time singularity is possible." 

Numerical simulations showed evidence of singular behavior when the ge- 
ometry of the level sets of the active scalar contain a hyperbolic saddle. There 
is a naturally associated notion of simple hyperbolic saddle breakdown. The 
main theorem we present in this paper shows that such breakdown cannot oc- 
cur in finite time. We also show that the angle of the saddle cannot close in 
finite time and it cannot be faster than a double exponential in time. Using 
the same techniques, we see that analogous results hold for incompressible 2D 
and 3D Euler. 

These results were announced in [9], but with a slight difference in the 
definition of a simple hyperbolic saddle. The definition given in Section 4 gen- 
eralizes the one given in the announcement. See also Constantin [4], discussed 
in Section 7, Remark 5 below. 

The whole work described in this paper is basically part of the author's 
thesis. I am particularly grateful to my thesis advisor Charles Fefferman for his 
attention, support, guidance and advice. I am indebted to D. Christodoulou 
and P. Constantin for helpful corrections and suggestions. I wish to thank 
A. Majda for suggesting the subject and E. Tabak for discussions and com- 
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ments on the first version of the paper. The author was supported by the Sloan 
Foundation during part of the writing of this paper. 



2. Quasi-geostrophic equation as a 2D model for 3D Euler 

The vorticity equation for incompressible 3D Euler is the following 
(1) (dt + u-V)u = (Vu)-oj, 

where uj is the vorticity and the velocity is divergence free. Using the Biot- 
Savart law we can recover the velocity from the vorticity by 

J_ r yxu(x + y,t) 
' ' br ./,,.. |y|S 

Beale, Kato and Majda [1] showed that a necessary and sufficient condition 
to have a singularity at time T is that 



L 



T 

\u{t)\\oodt = +00. 



'0 

Based on this result, Constantin, Fefferman and Majda [5] proved that if the 



direction field = is smooth in and near regions of high \co\ then 



blow-up does not occur. 

Another way to understand the problem is by constructing and studying 
models in lower dimensions. Constantin, Lax and Majda [6] studied a one- 
dimensional model for the 3D Euler vorticity equation, 

£ = H^. 

Here H is the Hilbert transform and the velocity is defined by u = J_ 00 co(y)dy. 
They show that there are solutions with breakdown. 

A 2D model of the quasi-geostrophic equation was studied in [7]. This 
equation represents a geophysical model [14], where the conjectured singular- 
ities describe the formation of sharp fronts between masses of hot and cold 
air. 

The 2D quasi-geostrophic active scalar equations are 

(2) (d t + u ■ v) e = o 

u = V^V where 6 = — (— A)^tp 

and the initial condition 0(x,0) = 9$. 

Here 9 = 9(x,t) with x G R 2 ,i G R + is a scalar, u is the velocity, and 
is the stream function. 
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In [7], [14], 9 is interpreted as a potential temperature. I am grateful to 
R. de la Llave for helpful comments, and for pointing out to me that in the 
quasi-geostrophic equation, 9 may also be regarded as a vorticity. 

In [7] it was shown that the level sets of 9 are analogous to the vortex lines 
for 3D Euler and that there is a geometric and analytic similarity between both 
equations. It was also shown there that results analogous to the ones mentioned 
above for 3D Euler also hold for the 2D active scalar equation: 

(i) By differentiating the equation we obtain an equation similar to (1) but 
in this case, instead of the vorticity, there is the gradient perpendicular to 
9: 

(d t + u- V) V ± 9 = (Vu) • V L 9. 
The stream function can be written 



x f 0(x + y) 



IB? 

thus 



i(x,t) = - 

JR 



V ± 9(x + y,t) 



dy, 



'i? 2 \y\ 

and obviously the velocity is divergence free. 

(ii) If 9 (x) € H k (~R 2 ) k > 3, then a necessary and sufficient condition for 
having a singularity at time T is that 

(3) / HV^Hoodt = +oo. 

Jo 

(iii) If the direction field = j^xfj is smooth in and near regions of high 
[V^l then blow-up does not occur. 

Another similarity is that in both cases the kinetic energy is conserved for 
all time. 



3. Numerical simulations 



It is not known if the quasi-geostrophic equation develops singularities in 
finite time. There have been several numerical attempts to find a candidate 
for initial data such that a strong singular behavior develops in finite time. In 
[7] the initial data 

9(x, 0) = sin(xi) sin(x2) + cos(a;2) 

were studied numerically and found to develop a strong front. An empirical 
asymptotic fit for maxjV^^I was obtained behaving like (T* — i) -1 - 66 with 
T* ^ 8.25. 
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2.5 5. 2.5 5. 

Figure 1. Numerical simulations 



Figure 1 represents the evolution of the initial data mentioned above at 
times t = 0,2, 4, 6. The lines are the level sets of 9. Is clear that for t = the 
level sets contain a hyperbolic saddle. As time evolves, Figure 1 shows that 
the saddle is closing very fast. Therefore there is a strong front formation. 
The author thanks P. Constantin, A. Majda and E. Tabak for permission to 
reproduce these pictures from their paper ([7]). 

Ohkitani and Yamada [13] gave another interpretation of this particular 
front. They suggested that the maximum gradient does not go to infinity in 
finite time, but rather goes to infinity at a double exponential rate in time. 
Recently, Constantin, Nie, and Schorghofer [8] made careful measurements of 
stretching rates and collected substantial numerical evidence against a singu- 
larity for this particular case. 

No other initial data that develop a front faster than the one mentioned 
above have been found. 
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4. Hyperbolic saddle scenario 

To rigorously and analytically approach the problem where the level sets 
of the active scalar field contain a hyperbolic saddle, we first need a precise 
definition of a hyperbolic saddle. 

In this section we define a naturally associated notion of simple hyperbolic 
blow-up scenario and we rule out a finite time singularity. 

Definition 1. A simple hyperbolic saddle in a neighborhood U of the origin 
is the set of curves p = const where 

p=( yi (3(t)+y 2 )( yi 6(t)-y 2 ), 

and there is a nonlinear time-dependent coordinate change 

yi = F 1 (x 1 ,x 2 ,t) 

2/2 = F 2 (xi,x 2 ,t) 

with p(t),5(t) e ^([0,^)), Ft e C 2 (U x [0,7;]), \0\, \6\ < C, /3(t) + S(t) > 0, 
|det§|| > c> whenever x£U,te [0, T*]. 

Remark 1. The saddle is allowed to rotate and dilate with respect to time. 
The center of the saddle can move in U with time. 

Remark 2. The angle of opening of the saddle is 7 ~ [3+5. The restrictions 
P(t) + 5(t) > and \(3\, \S\ < C are unimportant. See also [11]. 

Remark 3. The definition given in [9] for p was p = y\y 2 — cot a • y 2 2 where 
only one branch, y 2 = tan a ■ y\, was allowed to close and the other one moved 
smoothly. In the present paper both branches (y 2 = —f3(t)yi,y 2 = S(t)yi) are 
allowed to close at different rates. The results of [9] are therefore a particular 
case of the one presented here. 

The possible singularity in this scenario is due to j(t) becoming zero in 
finite time. The following theorem will show that this is not possible and j(t) 
can go to zero at most as a double exponential. 

Theorem 1. Let 6(x\,x 2 ,t) be a smooth solution of (2) defined for < 
t < T*, (x±,x 2 ) € R 2 ■ Assume for < t < T* that 9 is constant along the 
curves p = const defined in Definition 1. Assume also, for each fixed t, that 6 
is not constant on any disc in U. Then lim t ^T t 7(t) exists and is not 0. 

Corollary 1. Let 6 be as in Theorem 1, let £ = j^rfj 1 an d assume 
|V£| < C on (R 2 \ U) x [0,T*]. Then 6 continues to some solution of (2) on 
R 2 x [0, T* + e] for some e > 0. 
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Theorem 2. Let 0{x\,X2,t), (3, S, U and Fj be as in Theorem 1, but 
with T* = oo. Assume that the C 2 seminorms of Fj are bounded for all time 
t G [0,oo). Then 



log log — — I < (constant) • t + const 
7W 



for all t. 



Corollary 2. Let 9 be as in Theorem 2. Let £ = y^r§| • Assume |V£| < 
$(t) on (R 2 \U)x [0, oo). Then 



|V X ^| < e e J " 



on i? 2 x [0, oo). 



Remark 4. In view of Theorem 2, the most natural example is &(t) = e e 
Corollary 2 then shows that \V ± 9\ is bounded by a quadruple exponential. 



5. Stream function under a change of variables 



The purpose of this section is to obtain an expression for the stream 
function under a new set of variables (p,cr). The variable p = p(xi,X2,t) was 
defined in Definition 1, and a is defined as follows: 




y- Coordinate x- Coordinate 



Figure 2. Change of variables 

Let B(yi,y2,t) = be the bisector of the angle 7, p > and we denote 
4>(p) to be the intersection of B = with p as shown in Figure 2. We write 
<t>{x) for (f)(p(x)). 

We then define o~(x) to be the solution of 

(4) exp(aV^p)[4>(x)]=x. 

Note that, as a varies, exp(o"V^r/9)[0(x)] moves along the curve p = const. 
Thus a(x) is well defined for x G U and p > 0. 
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Let 0(xi,X2,t) be a smooth solution of (2) denned for < t < T*, 
(x\,X2) € R 2 . Assume for < t < T* that 9 is constant along the curves 
p = const where p = p(x±,X2,t). Hence we write 6(x, t) = 9(p, t) for a function 
9. Also, ip(x,t) = tp(p,a,t) and V L 9 can be expressed as 



where r = vCi+Cf an< ^ 

dx\ dp 

1 da dx2 
_ dx 2 _ dp 

2 da dxi ' 

By making the change of variables in equation (2), we obtain 

U .V6> = — (u-V ) = —( 9 ^ 9X2 I ^ 19x1 )■ 
<9p x 3p 9x2 9cr dx\ da ' 

90 d# /dp frp\ 
~di + 'd~p\di~^) ~ ' 



thus 



Taking into account that || is independent of cr and that 6 is not constant 
in a disk, we see that the derivative with respect to time of p along trajectories, 
i.e. ^ — is equal to a function (—Hi) that is also independent of a. 
Therefore 

< 5 > 

Furthermore, integrating with respect to a gives 

f 7 dp 

(6) iP(p,a,t)=H 1 (p,t)-a + j -£da + H 2 {p,t), 

where ^(p, t) = 4>(p, 0, i). 



6. Proofs 



The main idea in the proof of both theorems is to estimate the difference 
of the value of the stream function at a point p that lies in the branch of the 
saddle 1/2 = —j3(t)yi with the value of the stream function at q that lies in 
the other branch 7/2 = d(t)y\. Both p and q have the same y\ coordinate. 
We need two expressions of the stream function; one comes from the equality 
9 = — (— A)2ip and the other one comes from the change of variables done in 
Section 5. 
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Figure 3. Denning p, q, pi, qi, and r 

Lemma 1. Let 6 be a solution of equation (2), p and q be as defined above 
and ij) be given by ip = — (— A)~2#. Then 

\ip(p) - tp(q)\ < KiYf ■ log 7 1 

where K\ is a constant and \p — q\ ~ 7. 

To prove this lemma we use the fact that 9 is bounded for all time; this is 
because the derivative with respect to time along the Lagrangian trajectories 
is zero. We also recall that 0(x, t) belongs to I? for fixed t, with L 2 norm 
independent of t. 

i = - y,(<z) = / %)(— - Y^—My. 

\y-p\ \y-q\ 

Let us denote r = \p — q\. We split the integral /: 

I(x) = [ + [ + [ =h+I 2 + I 3 , 

J\h-p\< 2t ^2r<|j/-p|<fc Jk<\y-p\ 

where A; is a fixed number. 
We bound I\ by 

1 

?\<2r [ \y-p\ \y-Q 

1 

\y-p\ ' \y-Q 

< Ct. 



\h\<\\e\\ L ~- f It^-t^-Mv 

J\y-p\: 

<c- f (^^^ — -Wh 

J\y-p\<2r \y~ 
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We define s to be a point in the line between p and q. Then \y—p\ < 2\y—s\, 
and I2 can be estimated by 

1 

/ \v-s\ 
1 



\I 2 \<Ct- max s |V( , - 

J2T<\y-p\<_ 



< Ct ■ I max s - — - — pjdy 

hr<\y~p\<k \V ~ s \ 

< Ct ■ I logr|. 

To estimate I3 we recall that J R2 \9\ 2 dx is conserved for all time. It is easy 
to check that I/3I < C ■ r. Therefore |/| < r • | logr|. 

To state and prove the following lemmas we have to define 

q(yi,t) = {yi,5(t)y{) 

p(yi,t) = (yi, ~/3(t)yi). 

Lemma 2. With the same assumptions as in Theorem 1, let tp be given 
by expression (6) and (p, q) defined as before. Then 

where D = |det|^|. 

Proof. We evaluate ip at the points p\ = (p, a\) and q\ = (p, a 2 ) with 
G\ 7^ a 2 , as in Figure 3 

f 72 dp 

i>{<li) - i>{Pi) = Hi(p,t) ■ (cri - cr 2 ) + J —do. 

The next step is to take the limits p\ — > p and q\ — > (7. That means p — > 0. 

For the first term of the right-hand side of the equality we consider t £ 
[0, T*); hence the velocity is bounded as a function of x. By the construction 
of the set of variables (p, a) in Section 5, it follows that the cr-coordinates of 
p, q grow at most as log/? when p\ — > p and q± — ► (7. On the other hand .ffi is 
independent of a and is given by the following equation: 

For a fixed £ € [0, T*) we can bound Hi by 

|tfiOM)|<M-|V/9| + |^|, 

where 

dp _ dp_dm dp dy 2 

dxi dyi dxi dy 2 dxi 
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dp _ dp dyi dp dy 2 d(5 - (5) d{(58) 2 

dt ~ d Vl dt + dy 2 dt dt yi ' m+ dt Vl ' 

We take into account the definition of p = p(yi,y 2 ,t), y = F(xi, x 2 ,t) and 
Fi eC 2 (U x [0,T*]); then 

. dp 

\ — | < \y\ ■ (const) 

.dp. . . . . 
\ — \ < \y\ ■ (const). 



Thus 



\Hi{ P M < \y\- (const) 



(these constants may depend on t), and, from the definition of p, we also know 
that \y\ ~ pa when we approach the origin along the bisector B with p > 0. 
Therefore |i?i| is at most pa when p — > 0. This implies that 

lim Hi(p,t) ■ (oi - <j 2 ) = 0. 

For the second term we obtain something different. We define 

T = [(3/1,2/2) : P = const]. 
We have the following expressions for ^ by the change of variables 
% _ dyi dxi dyj_dx2 _ _dyi_dp_ dyi dp 



da dx\ da dx 2 da dx\ dx 2 dx 2 dx\ 

d yi ( dp dyi df^dy^A dm_ ( dp_dm_ dp dy 2 



dxi \dyidx 2 dy 2 dx 2 J dx 2 \dy\dxi dy 2 dx\ 



Then 



da dy 2 ' da dyi 

dyi da = 1 dy^da^ = 1 p 
<9cr dyi ' <9cr dy2 

Using these formulas for the integral on T we find that 

P ^l da = r ^^l d a + r ——da 

d(6-/3) f a2 , d(/35) r 2 2l 
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where 

dp dyi „ 



Jo D-M- 



f m 1 dy 2 



I: 



3 



-f 

Jo 

d(S - (3) /™ yl • y 2 



dt 7o d ■ M- 

d(p6) r y\ 2 



We obtain the following estimates for Ij(g) — 
Case /i. 



/3 <9yi 



» ft® j.r ^ a r ( ft® 



o{i). 



Case 12- 



h{q) ~ hip) = jT ( 
= 0( 7 ). 



Case I3. 



l D(p(y u t)) Diqiy\,t)) 



Case I4. 

Above we used the fact that 

^-(q(yi,t)) = -yiiP + 6), ^-(p(y 1 ,t)) = y 1 if3 + S) 

oy2 oy2 

^-(q(y~ 1 ,t)) = 6y 1 i(3 + 5), ^-(p( y ~ 1 ,t))={3y 1 iP + 5), 

oyi ayi 

and that the function is C 1 ^ x [0,T*]). 
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Now we can complete the proof of Lemma 2 by adding 7| to I| 

h+h -^t-J Dm,t)) dVl+ d[-J D(p(ui,t)) dVl - 

Lemma 3. Under the same assumptions as above and (q,r) as in Fig- 
ure 2, we then have 

s 2 = m - m = f t ■ £ Diitw) ^ 1 + E{xi ^ t] 

where E is bounded for all time. 

To prove the equality we again use the expression (6) for the stream func- 
tion, and evaluate it at 4>{p) and q\, where these points belong to the same level 
set p. Then we take the limit p — > 0, which means that <j>(p) — ► r and q\ — > q. 
The equality follows from using the same steps as in the proof of Lemma 2. 
The first two cases, I\ and I2, give us the function E(xi,X2,t), and the other 
two cases give us the other term of the right-hand side of the equality. 

Let us define the functions 

im= L mm"" 1 - 

Then Lemma 3 can be written like 



1C 

(7) —K(q) = S 2 -E(x 1 ,x 2 ,t), 

and from Lemma 2 



(8) S 1 = (- + -£)K(p) + -[K(q) - K(p)} + 0( 7 ). 

In (7) and (8), we know that S2 and E are bounded for all time and that 

c?7 dS d(3 
dt dt dt 

We also know that that there are two positive constants M and c such that 
M > K(q) > c> 0, M > K(p) > c> and 

K{p)-K{q)=0{ 1 ). 
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By combining (7) and (8) we obtain 

S 1 = ^K(p) + 0( 1 ). 

Now we use Lemma 1 to get the estimate 

(9) |^| < (const) | 7 - log 7 |, 

if 7 is less than a small constant. 

The proof of Theorems 1 and 2 follows directly from integrating (9). 



7. Stretching formula 

In the previous section we proved that if a solution 9 of equation (2) is 
constant along hyperbolas then the angle of the saddle 7 cannot close in finite 
time. The angle cannot close faster than a double exponential in time. 

In this section we want to prove Corollaries 1 and 2. Furthermore, if 
we assume that J* * &(s)ds is bounded for finite time then we can rule out a 
simple hyperbolic blow-up scenario by the necessary and sufficient condition 
(3) mentioned in Section 2. 

We use the formula 

(10) (^ + n ' v ) l v± 0| = «l v± 0|> 

where a is the level set stretching factor 

a= ^(Vn + V*u)£-£. 

Here, £ is the direction field of V^O. For (10), see Constantin [3] and Con- 
stantin, Majda and Tabak [7]. Consequently a can be represented by 

a(x) = P.V. / ((A • e(x)Mx + y) ■ ^x)))\V^9\(x + yA. 

Jb? \y\ \y\ 

In [7] they derive an upper bound for the magnitude of the stretching 
factor 

|«(x)| < C[G(x)\u{x)\ + (?G(x) + l)(G||e|| L =c + f- 2 ||^|| L2 )]. 
C, f are fixed constants and 

G{x) = sup| y |< f |V£(x + y)|. 
Now we derive an estimate for the velocity u(x): 



u(x) = I 
Je 



y^6{x + y) 



R2 1 2/ 1 3 



dy. 
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We consider r > and we define K(r) to be a smooth nonnegative function 
such that N(r) = 1 for < r < 1, and N(r) = for r > 2. Then we split u(cc) 
in two integrals: 

«(ar) =h + h = J N(— ) p + J (1 - N(— )) p dy. 

By integrating by parts, we can easily prove that 

< Ct ■su P \v ± e\ + c. 

For I2 we obtain 

J\y\>2r \y\ 2 

= 1 + / 

J2r<\y\<k J\y\>k 

where k is fixed. 

By performing the first integral and for the second one using the conser- 
vation of J R2 \6\ 2 dx, 

\h\ < C\ logr| + const. 

Therefore 

\u(x)\ < CisuplV^^lr + C 2 \ log -| + C 3 . 

r 

Take r = sup |y± e | ■ Then u is bounded by logsuplV- 1 ^, assuming that 
suplV- 1 ^ > e. 

Under the assumptions of Theorem 1, G{x) is bounded by a double ex- 
ponential for \x\ > c. Therefore we can estimate a on V = U \ {\x\ < c} 
by 



(11) \a\ < logHV^^IIiooe 

Proposition 1. Assume the existence of a function f = f(p,t) inde 

dt 



pendent of a. Then the material derivative (= D t = jL + u ■ V) of f is also 



independent of a. 

To prove the proposition we have to compute the material derivative: 
, ( 9 „\ r/ , df df (dp dip 

From equation (5) and the fact that the functions |^ are independent of 
<r, we can conclude that the material derivative of / is also independent of a. 
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In particular the function F 
material derivative of F as follows 

1 



(12) 



|Vp| ' |Vp| 



|y^l is independent of <r. We write the 
A(|V^|) + |V^|A( ' 



|Vp|" 



= [a + \Vp\D t {- 



1 



|Vp|V |v P | 



We can estimate \V p\D t {^-^) on F by 
(13) | V p||A(^j)| = ^|A(|Vp|)| 



■|Vp|- 



|Vp| 

1 



9* 



,d<y, 



+ u- V ) |Vp| 
.dp. 



d/3\ 



^Cil^l + Csl^l + ClogHV^^IlLooe 
at at 

Lemma 3 shows that |^| < const. The same estimate can be obtained, |^ 
< const, by substituting the point q with p in Lemma 3. 
Combining (11) and (13) we obtain 

1 



« + |Vp|A(^) 



< e e log | |V 9 1| L°o + const 



on V when HV^lc, > e e \ Then by (12) 

(14) \D t F\ <e et log||V^|| L oo F 

on y when HV^^Hloo > e e \ However, because F only depends on p and t and 
by the proposition it follows that inequality (14) only depends on p and t. 

If |x| < c, then find x € V with p(x) = p(x). Inequality (14) holds for x 
and therefore it holds for x. Thus, (14) holds, not just on V, but on U. 

We want to estimate the material derivative of (V^l on U. By (12) we 
have the following equality 

(15) A(|V^|) = |Vp| • D t F + F • A(|Vp|). 

The estimate for the first term of the right-hand side of equality (15) is obtained 
in (14). For the second term we know that on V, 

(16) F < ||V~ l 0||loo. 

Because F is independent of a, inequality (16) holds on U. We compute the 
material derivative of |Vp| and estimate it on U as follows: 

d 



|A(|Vp|)| 



'dt 



|Vp| + u-V(|Vp|)| 



^^l + ^lfl + l-vdvpl)! 

. „ , do . . dp . ._i .„ P t 

< C 1 \ — \ + C 2 \-£\ + log\\V J -9\\ L o e e . 
dt dt 
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We recall again Lemma 3 and get the following estimate on U: 

lAdV^DI < CHV^Hloo log HV^IlLcoe 6 '. 

If x e R 2 \ U then |V£| < $(£). By using the upper bound of a we obtain 

|A(|V^|)| <C||V^|| L oc log||V^|| L oc(d>(t) + e et ) 

on x € R 2 \ U. 
Therefore 

(17) I^HV^HlocI < CHV^IIlcc log ||V^|| L oc(e et + $(t)) 

whenever ||V j -^||loo > e e \ 

At time t = 0, ||V~ l 0||loo < oo. We finish the proofs of Corollaries 1 and 
2 by integrating (17). 

Remark 5. In [4], Constantin shows that under the hypothesis of a "proper 
nondegenerate-self-similar Ansatz" the ||V- l 6 , ||lcx) grows at most as a double 
exponential. 

Remark 6. Suppose 9 = 9(H(x,t),t) is constant along ellipses. We allow 
these ellipses to close according to 

U(x,t) = a{t) ■ x\ + b{t) ■ x\. 

Whenever this happens in the numerical simulations, the norm of the gradient 
of 9 does not grow fast. We can show that in this case the (V -1 ^! is bounded 
by a double exponential in time. 

8. Similar results for Euler 

Two dimensional case: Majda and Tabak [12] ran numerical simulations 
with the same initial data for the 2D Euler vorticity equation and compared 
them with the simulations in [7]. The norm of the gradient perpendicular of 
the vorticity grew much slower for 2D Euler than the \V ± 9\ above. 

The representation of 2D Euler equation in vorticity form is 

(d t + u ■ V) co = 
u = V ±r ip where u = Aip. 

The two active scalars 9 and u> are similar, but they differ in the characteriza- 
tion of the stream function. 

Using the same scheme as in Section 4, we assume that oj is constant along 
hyperbolas and that 7 is the angle of the saddle. We can show 

I log 7(i) I < (constant) • t, 
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which means 7 can go to zero at most as an exponential. The proof is identical 
to the one in Section 6, but in Lemma 1 tp is defined by 

V> = — / u(x + y) log \y\dy 
27T Jr? 

and it follows that 

hKp)-#?)| <K\j\. 

Three-dimensional case: As we explained in Section 2, the quasi-geostro- 
phic equation is a two-dimensional model for the 3D incompressible Euler 
equation. The techniques used in this paper give analogous results for 3D 
Euler [10]: 

Theorem 3. Let u(x, t) be a smooth solution of the 3D Euler incom- 
pressible equation defined for < t < T*, x £ R 3 with 

\u\ dp dp 
uj = curl (u) = — (- ^— , j—, . 

r 0x2 ox 1 

Here, r = |Vp|, u is bounded up to t = T* and p is defined as in Theorem 1 
with the same nonlinear time dependent coordinate change and the same as- 
sumptions. If uj is not zero in a disc in U, then lim^y* ^{t) exists, and is 
not 0. 

Corollary 3. Let u be as in Theorem 3. Let C = o- Assume |V£| < C 
on (i? 3 \ U) x [0, T*]. T/ien u continues to some solution of the 3D Euler 
incompressible equation on R 3 x [0, T* + e] for some e > 0. 

A similar estimate is obtained for the angle of the saddle 7. Note that 
vortex stretching may take place in the setting of Theorem 3. 

Princeton University, Princeton, NJ 
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